A Novel Mathematical Model That Predicts the Protection Time of SARS-CoV-2 Antibodies.

Infectious diseases such as SARS-CoV-2 pose a considerable threat to public health. Constructing a reliable mathematical model helps us quantitatively explain the kinetic characteristics of antibody-virus interactions. A novel and robust model is developed to integrate antibody dynamics with virus dynamics based on a comprehensive understanding of immunology principles. This model explicitly formulizes the pernicious effect of the antibody, together with a positive feedback stimulation of the virus-antibody complex on the antibody regeneration. Besides providing quantitative insights into antibody and virus dynamics, it demonstrates good adaptivity in recapturing the virus-antibody interaction. It is proposed that the environmental antigenic substances help maintain the memory cell level and the corresponding neutralizing antibodies secreted by those memory cells. A broader application is also visualized in predicting the antibody protection time caused by a natural infection. Suitable binding antibodies and the presence of massive environmental antigenic substances would prolong the protection time against breakthrough infection. The model also displays excellent fitness and provides good explanations for antibody selection, antibody interference, and self-reinfection. It helps elucidate how our immune system efficiently develops neutralizing antibodies with good binding kinetics. It provides a reasonable explanation for the lower SARS-CoV-2 mortality in the population that was vaccinated with other vaccines. It is inferred that the best strategy for prolonging the vaccine protection time is not repeated inoculation but a directed induction of fast-binding antibodies. Eventually, this model will inform the future construction of an optimal mathematical model and help us fight against those infectious diseases.

More or less deadly? A mathematical model that predicts SARS-CoV-2 evolutionary direction.

SARS-CoV-2 has caused tremendous deaths globally. It is of great value to predict the evolutionary direction of SARS-CoV-2. In this paper, we proposed a novel mathematical model that could predict the evolutionary trend of SARS-CoV-2. We focus on the mutational effects on viral assembly capacity. A robust coarse-grained mathematical model is constructed to simulate the virus dynamics in the host body. Both virulence and transmissibility can be quantified in this model. A delicate equilibrium point that optimizes the transmissibility can be numerically obtained. Based on this model, the virulence of SARS-CoV-2 might further decrease, accompanied by an enhancement of transmissibility. However, this trend is not continuous; its virulence will not disappear but remains at a relatively stable range. A virus assembly model which simulates the virus packing process is also proposed. It can be explained why a few mutations would lead to a significant divergence in clinical performance, both in the overall particle formation quantity and virulence. This research provides a novel mathematical attempt to elucidate the evolutionary driving force in RNA virus evolution.

Statistical analysis supports UTR (untranslated region) deletion theory in SARS-CoV-2.

It was noticed that the mortality rate of SARS-CoV-2 infection experienced a significant declination in the early stage of the epidemic. We suspect that the sharp deterioration of virus toxicity is related to the deletion of the untranslated region (UTR) of the virus genome. It was found that the genome length of SARS-CoV-2 engaged a significant truncation due to UTR deletion after a mega-sequence analysis. Sequence similarity analysis further indicated that short UTR strains originated from its long UTR ancestors after an irreversible deletion. A good correlation was discovered between genome length and mortality, which demonstrated that the deletion of the virus UTR significantly affected the toxicity of the virus. This correlation was further confirmed in a significance analysis of the genetic influence on the clinical outcomes. The viral genome length of hospitalized patients was significantly more extensive than that of asymptomatic patients. In contrast, the viral genome length of asymptomatic was considerably longer than that of ordinary patients with symptoms. A genome-level mutation scanning was performed to systematically evaluate the influence of mutations at each position on virulence. The results indicated that UTR deletion was the primary driving force in alternating virus virulence in the early evolution. In the end, we proposed a mathematical model to explain why this UTR deletion was not continuous.

A Continuous Markov-Chain Model for the Simulation of COVID-19 Epidemic Dynamics.

To address the urgent need to accurately predict the spreading trend of the COVID-19 epidemic, a continuous Markov-chain model was, for the first time, developed in this work to predict the spread of COVID-19 infection. A probability matrix of infection was first developed in this model based upon the contact frequency of individuals within the population, the individual's characteristics, and other factors that can effectively reflect the epidemic's temporal and spatial variation characteristics. The Markov-chain model was then extended to incorporate both the mutation effect of COVID-19 and the decaying effect of antibodies. The developed comprehensive Markov-chain model that integrates the aforementioned factors was finally tested by real data to predict the trend of the COVID-19 epidemic. The result shows that our model can effectively avoid the prediction dilemma that may exist with traditional ordinary differential equations model, such as the susceptible-infectious-recovered (SIR) model. Meanwhile, it can forecast the epidemic distribution and predict the epidemic hotspots geographically at different times. It is also demonstrated in our result that the influence of the population's spatial and geographic distribution in a herd infection event is needed in the model for a better prediction of the epidemic trend. At the same time, our result indicates that no simple derivative relationship exists between the threshold of herd immunity and the virus basic reproduction number R0. The threshold of herd immunity achieved through natural immunity is significantly higher than 1 - 1/R0. These not only explain the theoretical misconceptions of herd immunity thresholds in herd immunity theory but also provide a guidance for predicting the optimal vaccination coverage. In addition, our model can predict the temporal and spatial distribution of infections in different epidemic waves. It is implied from our model that it is challenging to eradicate COVID-19 in the short term for a large population size and a wide spatial distribution. It is predicted that COVID-19 is likely to coexist with humans for a long time and that it will exhibit multipoint epidemic effects at a later stage. The statistical evidence is consistent with our prediction and strongly supports our modeling results.